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We present results for some of the light baryon masses and their excited states in 2+1 flavour 
domain wall QCD. We considered several lattice spacings, with the DBW2 and Iwasaki gauge 
actions and different sea quark masses on a volume of 16^ x 32 and a fifth dimension of size 8. 
All data were generated on the QCDOC machines. Despite large residual massses and a limited 
number of sea quark mass values with which to perform chiral extrapolations, our results are 
in reasonable agreement with experiment and scale within errors. Finite size effects on most 
ensembles appear to be smaU. 
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1. Introduction 

The calculation of many quantities from lattice QCD with quenched, or heavy dynamical, 
quarks has produced results which, although in qualitative agreement with experiment, have un- 
controlled sources of error. To progress beyond this, calculations with light dynamical quarks are 
necessary. Using fermions which satisfy the Ginsparg-Wilson (GW) relation is advantageous, as 
the theory is well defined within broad parameter ranges. It has lattice versions of the symmetries 
of continuum QCD, and the correct flavour content. Mixing of operators with different chirality is 
suppressed, renormalisation is simplified, and a continuum-like chiral perturbation theory can be 
used for the extrapolation of quantities to the chiral limit. Domain wall fermions (DWF) satisfy 
these requirements. The QCDOC machine and the RHMC algorithm have allowed for the first time 
calculations using DWF with 2 + 1 flavours of dynamical quarks. 

Baryon physics is an area rich in phenomenology. Unanswered questions as to the nature of 
some of the excited states, the decay of the proton predicted by some Grand Unified and super sym- 
metric models, and the determination of matrix elements related to the structure functions and the 
neutron electric dipole moment can, in principle, be determined by lattice QCD calculations. These 
calculations are very computationally challenging, as they need both light sea quarks and large vol- 
umes, as well as fermions with the correct symmetry properties. In this work we report on a study 
of the lowest lying states in the baryon spectrum, {N ,^,Q.,N*}, on ensembles of configurations 
produced by QCDOC. It is important for any calculation claiming to be "full QCD" to be able to 
reproduce this spectrum. These initial ensembles were produced primarily for a search of parame- 
ter space to guide larger production runs. As such, they have a small volume, too small certainly for 
excited states such as the N*, but this can be used, eventually in combination with the productions 
runs on larger volumes to try and estimate the size of the finite volume effects on the remaining 
spectrum. 

The DWF action is given in |[T]] with the Pauli-Villars field in for the dynamical simulation. 
The gauge fields were generated with renormalisation group (RG) improved actions, as follows: 



with either ci = -1.4069 for the DBW2 action [|]] or ci = -0.331 for the Iwasaki action [@]. It 
has been noted in the quenched DWF calculations that these actions reduce the chiral symmetry 
breaking resulting from the finite fifth dimension. The mechanism and size of chiral symmetry 
breaking on these ensembles is studied in detail in 

In the DWF formalism, the 4D quark fields are constructed from left and right projections of 
the 5D fermion fields on the boundaries. With a finite fifth dimension, there is still an overlap 
between these left- and right-handed fields, which manifests itself as an additive mass renormalisa- 
tion, known as the residual mass. This can be determined directly from the axial Ward-Takahashi 
identity 



Sg = -^ (1-8)ci £pW^v + ci 



(1.1) 



d^A^{x) = 2amfP{x) +2Ji{x) 
Ri 2{amf + am^gs)P{x) 



(1.2) 
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j3 ami/am^ r^/a a HCeV) L(fm) Lmp mpjmy #traj 



0.72 
0.72 


0.04/0.04 
0.02/0.04 


4.3(lj 


1.7(1) 


1 o 

i.y 




7.7(1) 
6.0(1) 


0.692(5) 
0.589(3) 


3400 
6000 


0.764 
0.764 


0.04/0.04 
0.02/0.04 


5.1(2) 


2.0(1) 


1.6 


;i) 


6.7(1) 
5.1(1) 


0.699(4) 
0.619(4) 


5750 
3000 


2.13 
2.13 


0.04/0.04 
0.02/0.04 


4.6(2) 


1.8(1) 


1.8 


;i) 


7.5(1) 
5.8(1) 


0.700(8) 
0.615(5) 


3600 
3600 


2.2 
2.2 


0.04/0.04 
0.02/0.04 


5.3(1) 


2.1(1) 


1.5 


;i) 


6.8(1) 
5.1(1) 


0.726(2) 
0.667(8) 


4500 
3200 



Table 1: Properties of the ensembles used in this study. The value of ro/a was determined in [^. The lattice 
spacing, and thus the volume are set by choosing = 0.5fm. Trajectories are of length 0.5 

where 75 (x) is the point-split current constructed from fields at the mid-point of the fifth dimension 
and P{x) the pseudoscalar density. 

2. Details of the calculation 

The ensembles were all generated on QCDOC machines, using the exact RHMC algorithm [^, 
|8]], with a volume of 16^ x 32 and a fifth dimension of size Ls = 8 and are detailed in Table [l[ The 
integrated autocorrelation times are estimated to be G (50) for mesonic correlators. To maximise 
statistics, the correlation functions were measured every 5 trajectories, and then binned so that the 
separation between independent measurements is 100 trajectories. The correlators were computed 
from up to four sources on different time-planes with as many as three different smeared sources. 

The masses of the pseudoscalar and vector mesons were determined by simultaneous fits to 
the local and smeared correlators with two exponentials, the ground state and the excited state. 
Similarly /Wres was determined by a simultaneous fit of a plateau to a ratio of correlators, with both 
smeared and local sources. 

The standard baryon interpolating operator is given by 

= £ijk [Wi{x)CT\ifj{x)] Yk{x) (2.1) 

For the / = |} baryons, F = {75,757*:}- Another operator, which projects onto the negative 
parity I = j state, with F = 1 was also used. For baryon correlators in a finite box with periodic 
boundary conditions, the backward propagating state is the negative parity partner, that is 

Cb(0 =A+e-"'^' +A-e-"'-''^-'^ (2.2) 

For the / = 5 baryon, the masses of the positive and negative parity states were determined by 
a simultaneous fit to equation ( p!^ ) using the standard operator, and a single exponential to the 
negative parity correlator. This is shown in Figure |I](a). Typically this was computed for the 
smeared correlator only, as the local correlator had a poor signal. 

The mass of the I = j baryon was determined from a single exponential fit to the smeared 
correlator, as the signal for the excited, negative parity partner was poor, as might be expected on 
relatively small numbers of configurations. 
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Figure 1: (a) The effective mass of / = ^ baryon coiTelator on tlie DBW2 j3 = 0.72, ami / am^ — 0.02 /0.04 
ensemble. The red symbol shows the data for the positive parity state, the closed blue squares are the time- 
reversed backward propagating negative parity state data and the open blue squares are the negative parity 
data from the second correlator (b) The Edinburgh plot for the different ensembles. 



3. Preliminary Results 

Shown in Figure ^b) is the Edinburgh plot. It is reassuring that, even at relatively coarse 
lattice spacing, with a small fifth dimension and consequently moderate chiral symmetry breaking, 
the data follows the phenomenological curve very well. The only exception is the lightest Iwasaki 
j8 = 2.2 datum. Naively, one might expect this to be a finite size effect, especially, examining 
Table 0, given that the lattice spacing and thus the box size are rather small. However, the value 
of Lmp is not significantly smaller than the other data sets and, critically, Lmp > 4, suggesting that 
the box is big enough as measured by the pseudoscalar meson. 

For this particular ensemble the signal for the vector meson mass is not good. In particular, the 
effective mass plot has a poor plateau, and a stable fit can only be achieved for the lowest region. 
The net effect is for a rather low vector meson mass. It is probable that this is due simply to low 
statistics. A low estimate of the vector mass would cause this datum in the Edinburgh plot to be 
shifted up and to the right, which could be mistaken for a finite volume effect increasing the mass 
of the nucleon. 

The quark mass was defined as 

aiTiq = amf + am^esifnf) (3.1) 

with the chiral limit at airiq = 0. With two sea quark masses, only a crude chiral extrapolation could 
be attempted, i.e. drawing a straight line through the two data points. The strange quark mass, (for 
the n baryon) was set from the kaon mass. 

Shown in Figure ^ is the dependence of the baryon spectrum, in dimensionless units, on the 
lattice spacing. A continuum extrapolation cannot be attempted with these ensembles. However, 
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Figure 2: Scaling of the baryon spectrum with lattice spacing squared. The symbols denote the following: 
Closed circles A^, open circles A^*, squares A, diamonds H. Black symbols denote experiment, red DBW2 
ensembles, blue Iwasaki ensembles. The value of ro = 0.5fm was chosen, to give an indication of the 
experimental spectrum in these units. The A^* measured on the j3 =2.13 ensemble (the furthest right open 
blue circle) has been offset to the left for clarity. 



the data for the ground states, {A/^, A,n}, shows reasonable scaling, albeit with large errors. These 
large uncertainties are due to the crude nature of the chiral extrapolation. The determination of ro is 
well defined for each ensemble, and thus it is a good quantity to use to examine scaling behaviour. 
As the value of ro in the continuum is unknown, rather than setting the absolute scale with ro, 
a better strategy is to predict dimensionless ratios of physical quantities in the continuum, while 
using ro to just examine the scaling behaviour. 

The negative parity partner of the nucleon, the A'^* is expected to become degenerate with the 
nucleon in a small enough box. This effect can be clearly seen from Figure ^; the A^* mass drops 
dramatically as the volume is reduced. This suggests that finite size effects may also be beginning 
to affect the ground states for the ensembles at finest lattice spacing. These finite size effects would 
tend increase the mass of the ground states. The slight upward tendancy in the scaling plot as 
the lattice spacing decreases is consistent with the finite size effects spoiling otherwise very good 
scaling, or a small scaling violation for the nucleon mass. 



4. Conclusions 

We have determined the spectrum of the lowest lying baryon states for several ensembles 
with two different gauge actions. The QCDOC machines and the RHMC algorithm have made 
2 + 1 flavours of DWF ensembles possible for the first time. With limited statistics, and only two 
different sea quark masses, we qualitatively reproduce the experimental spectrum. There is limited 
evidence that finite size effects may be influencing the ground state baryons on the smaller volumes. 
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Both the Edinburgh plot, and the scaling analysis suggest that a programme of baryon physics on 
larger volumes and at lighter quark mass will yield very interesting results. 
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